Comparative analysis of mitochondrial genomes of Schisandra repanda and Kadsura japonica.

The family Schisandraceae is a basal angiosperm plant group distributed in East and Southeast Asia and includes many medicinal plant species such as Schisandra chinensis. In this study, mitochondrial genomes (mitogenomes) of two species, Schisandra repanda and Kadsura japonica, in the family were characterized through de novo assembly using sequencing data obtained with Oxford Nanopore and Illumina sequencing technologies. The mitogenomes of S. repanda were assembled into one circular contig (571,107 bp) and four linear contigs (10,898-607,430 bp), with a total of 60 genes: 38 protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The mitogenomes of K. japonica were assembled into five circular contigs (211,474-973,503 bp) and three linear contigs (8,010-72,712 bp), with a total of 66 genes: 44 PCGs, 19 tRNA genes, and 3 rRNA genes. The mitogenomes of the two species had complex structural features with high repeat numbers and chloroplast-derived sequences, as observed in other plant mitogenomes. Phylogenetic analysis based on PCGs revealed the taxonomical relationships of S. repanda and K. japonica with other species from Schisandraceae. Finally, molecular markers were developed to distinguish between S. repanda, K. japonica, and S. chinensis on the basis of InDel polymorphisms present in the mitogenomes. The mitogenomes of S. repanda and K. japonica will be valuable resources for molecular and taxonomic studies of plant species that belong to the family Schisandraceae.

Long-read transcriptome sequencing provides insight into lignan biosynthesis during fruit development in Schisandra chinensis.

BACKGROUND: Schisandra chinensis, an ancient member of the most basal angiosperm lineage which is known as the ANITA, is a fruit-bearing vine with the pharmacological effects of a multidrug system, such as antioxidant, anti-inflammatory, cardioprotective, neuroprotective, anti-osteoporosis effects. Its major bioactive compound is represented by lignans such as schisandrin. Molecular characterization of lignan biosynthesis in S. chinensis is of great importance for improving the production of this class of active compound. However, the biosynthetic mechanism of schisandrin remains largely unknown. RESULTS: To understand the potential key catalytic steps and their regulation of schisandrin biosynthesis, we generated genome-wide transcriptome data from three different tissues of S. chinensis cultivar Cheongsoon, including leaf, root, and fruit, via long- and short-read sequencing technologies. A total of 132,856 assembled transcripts were generated with an average length of 1.9 kb and high assembly completeness. Overall, our data presented effective, accurate gene annotation in the prediction of functional pathways. In particular, the annotation revealed the abundance of transcripts related to phenylpropanoid biosynthesis. Remarkably, transcriptome profiling during fruit development of S. chinensis cultivar Cheongsoon revealed that the phenylpropanoid biosynthetic pathway, specific to coniferyl alcohol biosynthesis, showed a tendency to be upregulated at the postfruit development stage. Further the analysis also revealed that the pathway forms a transcriptional network with fruit ripening-related genes, especially the ABA signaling-related pathway. Finally, candidate unigenes homologous to isoeugenol synthase 1 (IGS1) and dirigent-like protein (DIR), which are subsequently activated by phenylpropanoid biosynthesis and thus catalyze key upstream steps in schisandrin biosynthesis, were identified. Their expression was increased at the postfruit development stage, suggesting that they may be involved in the regulation of schisandrin biosynthesis in S. chinensis. CONCLUSIONS: Our results provide new insights into the production and accumulation of schisandrin in S. chinensis berries and will be utilized as a valuable transcriptomic resource for improving the schisandrin content.

Systemic Expression of Genes Involved in the Plant Defense Response Induced by Wounding in Senna tora.

Wounds in tissues provide a pathway of entry for pathogenic fungi and bacteria in plants. Plants respond to wounding by regulating the expression of genes involved in their defense mechanisms. To analyze this response, we investigated the defense-related genes induced by wounding in the leaves of Senna tora using RNA sequencing. The genes involved in jasmonate and ethylene biosynthesis were strongly induced by wounding, as were a large number of genes encoding transcription factors such as ERFs, WRKYs, MYBs, bHLHs, and NACs. Wounding induced the expression of genes encoding pathogenesis-related (PR) proteins, such as PR-1, chitinase, thaumatin-like protein, cysteine proteinase inhibitor, PR-10, and plant defensin. Furthermore, wounding led to the induction of genes involved in flavonoid biosynthesis and the accumulation of kaempferol and quercetin in S. tora leaves. All these genes were expressed systemically in leaves distant from the wound site. These results demonstrate that mechanical wounding can lead to a systemic defense response in the Caesalpinioideae, a subfamily of the Leguminosae. In addition, a co-expression analysis of genes induced by wounding provides important information about the interactions between genes involved in plant defense responses.

Machine learning, transcriptome, and genotyping chip analyses provide insights into SNP markers identifying flower color in Platycodon grandiflorus.

Bellflower is an edible ornamental gardening plant in Asia. For predicting the flower color in bellflower plants, a transcriptome-wide approach based on machine learning, transcriptome, and genotyping chip analyses was used to identify SNP markers. Six machine learning methods were deployed to explore the classification potential of the selected SNPs as features in two datasets, namely training (60 RNA-Seq samples) and validation (480 Fluidigm chip samples). SNP selection was performed in sequential order. Firstly, 96 SNPs were selected from the transcriptome-wide SNPs using the principal compound analysis (PCA). Then, 9 among 96 SNPs were later identified using the Random forest based feature selection method from the Fluidigm chip dataset. Among six machines, the random forest (RF) model produced higher classification performance than the other models. The 9 SNP marker candidates selected for classifying the flower color classification were verified using the genomic DNA PCR with Sanger sequencing. Our results suggest that this methodology could be used for future selection of breeding traits even though the plant accessions are highly heterogeneous.

Chromosomal Mapping of Tandem Repeats Revealed Massive Chromosomal Rearrangements and Insights Into Senna tora Dysploidy.

Tandem repeats can occupy a large portion of plant genomes and can either cause or result from chromosomal rearrangements, which are important drivers of dysploidy-mediated karyotype evolution and speciation. To understand the contribution of tandem repeats in shaping the extant Senna tora dysploid karyotype, we analyzed the composition and abundance of tandem repeats in the S. tora genome and compared the chromosomal distribution of these repeats between S. tora and a closely related euploid, Senna occidentalis. Using a read clustering algorithm, we identified the major S. tora tandem repeats and visualized their chromosomal distribution by fluorescence in situ hybridization. We identified eight independent repeats covering ~85 Mb or ~12% of the S. tora genome. The unit lengths and copy numbers had ranges of 7-5,833 bp and 325-2.89 × 106, respectively. Three short duplicated sequences were found in the 45S rDNA intergenic spacer, one of which was also detected at an extra-NOR locus. The canonical plant telomeric repeat (TTTAGGG)n was also detected as very intense signals in numerous pericentromeric and interstitial loci. StoTR05_180, which showed subtelomeric distribution in Senna occidentalis, was predominantly pericentromeric in S. tora. The unusual chromosomal distribution of tandem repeats in S. tora not only enabled easy identification of individual chromosomes but also revealed the massive chromosomal rearrangements that have likely played important roles in shaping its dysploid karyotype.

Transcriptome Changes Reveal the Molecular Mechanisms of Humic Acid-Induced Salt Stress Tolerance in Arabidopsis.

Humic acid (HA) is a principal component of humic substances, which make up the complex organic matter that broadly exists in soil environments. HA promotes plant development as well as stress tolerance, however the precise molecular mechanism for these is little known. Here we conducted transcriptome analysis to elucidate the molecular mechanisms by which HA enhances salt stress tolerance. Gene Ontology Enrichment Analysis pointed to the involvement of diverse abiotic stress-related genes encoding HEAT-SHOCK PROTEINs and redox proteins, which were up-regulated by HA regardless of salt stress. Genes related to biotic stress and secondary metabolic process were mainly down-regulated by HA. In addition, HA up-regulated genes encoding transcription factors (TFs) involved in plant development as well as abiotic stress tolerance, and down-regulated TF genes involved in secondary metabolic processes. Our transcriptome information provided here provides molecular evidences and improves our understanding of how HA confers tolerance to salinity stress in plants.

De novo Transcriptome Assembly of Senna occidentalis Sheds Light on the Anthraquinone Biosynthesis Pathway.

Senna occidentalis is an annual leguminous herb that is rich in anthraquinones, which have various pharmacological activities. However, little is known about the genetics of S. occidentalis, particularly its anthraquinone biosynthesis pathway. To broaden our understanding of the key genes and regulatory mechanisms involved in the anthraquinone biosynthesis pathway, we used short RNA sequencing (RNA-Seq) and long-read isoform sequencing (Iso-Seq) to perform a spatial and temporal transcriptomic analysis of S. occidentalis. This generated 121,592 RNA-Seq unigenes and 38,440 Iso-Seq unigenes. Comprehensive functional annotation and classification of these datasets using public databases identified unigene sequences related to major secondary metabolite biosynthesis pathways and critical transcription factor families (bHLH, WRKY, MYB, and bZIP). A tissue-specific differential expression analysis of S. occidentalis and measurement of the amount of anthraquinones revealed that anthraquinone accumulation was related to the gene expression levels in the different tissues. In addition, the amounts and types of anthraquinones produced differ between S. occidentalis and S. tora. In conclusion, these results provide a broader understanding of the anthraquinone metabolic pathway in S. occidentalis.

Genome-enabled discovery of anthraquinone biosynthesis in Senna tora.

Senna tora is a widely used medicinal plant. Its health benefits have been attributed to the large quantity of anthraquinones, but how they are made in plants remains a mystery. To identify the genes responsible for plant anthraquinone biosynthesis, we reveal the genome sequence of S. tora at the chromosome level with 526 Mb (96%) assembled into 13 chromosomes. Comparison among related plant species shows that a chalcone synthase-like (CHS-L) gene family has lineage-specifically and rapidly expanded in S. tora. Combining genomics, transcriptomics, metabolomics, and biochemistry, we identify a CHS-L gene contributing to the biosynthesis of anthraquinones. The S. tora reference genome will accelerate the discovery of biologically active anthraquinone biosynthesis pathways in medicinal plants.

Humic acid enhances heat stress tolerance via transcriptional activation of Heat-Shock Proteins in Arabidopsis.

Humic acid (HA) is composed of a complex supramolecular association and is produced by humification of organic matters in soil environments. HA not only improves soil fertility, but also stimulates plant growth. Although numerous bioactivities of HA have been reported, the molecular evidences have not yet been elucidated. Here, we performed transcriptomic analysis to identify the HA-prompted molecular mechanisms in Arabidopsis. Gene ontology enrichment analysis revealed that HA up-regulates diverse genes involved in the response to stress, especially to heat. Heat stress causes dramatic induction in unique gene families such as Heat-Shock Protein (HSP) coding genes including HSP101, HSP81.1, HSP26.5, HSP23.6, and HSP17.6A. HSPs mainly function as molecular chaperones to protect against thermal denaturation of substrates and facilitate refolding of denatured substrates. Interestingly, wild-type plants grown in HA were heat-tolerant compared to those grown in the absence of HA, whereas Arabidopsis HSP101 null mutant (hot1) was insensitive to HA. We also validated that HA accelerates the transcriptional expression of HSPs. Overall, these results suggest that HSP101 is a molecular target of HA promoting heat-stress tolerance in Arabidopsis. Our transcriptome information contributes to understanding the acquired genetic and agronomic traits by HA conferring tolerance to environmental stresses in plants.

Whole-genome, transcriptome, and methylome analyses provide insights into the evolution of platycoside biosynthesis in Platycodon grandiflorus, a medicinal plant.

Triterpenoid saponins (TSs) are common plant defense phytochemicals with potential pharmaceutical properties. Platycodon grandiflorus (Campanulaceae) has been traditionally used to treat bronchitis and asthma in East Asia. The oleanane-type TSs, platycosides, are a major component of the P. grandiflorus root extract. Recent studies show that platycosides exhibit anti-inflammatory, antiobesity, anticancer, antiviral, and antiallergy properties. However, the evolutionary history of platycoside biosynthesis genes remains unknown. In this study, we sequenced the genome of P. grandiflorus and investigated the genes involved in platycoside biosynthesis. The draft genome of P. grandiflorus is 680.1 Mb long and contains 40,017 protein-coding genes. Genomic analysis revealed that the CYP716 family genes play a major role in platycoside oxidation. The CYP716 gene family of P. grandiflorus was much larger than that of other Asterid species. Orthologous gene annotation also revealed the expansion of ß-amyrin synthases (bASs) in P. grandiflorus, which was confirmed by tissue-specific gene expression. In these expanded gene families, we identified key genes showing preferential expression in roots and association with platycoside biosynthesis. In addition, whole-genome bisulfite sequencing showed that CYP716 and bAS genes are hypomethylated in P. grandiflorus, suggesting that epigenetic modification of these two gene families affects platycoside biosynthesis. Thus whole-genome, transcriptome, and methylome data of P. grandiflorus provide novel insights into the regulation of platycoside biosynthesis by CYP716 and bAS gene families.

Genome Assembly and Annotation of Soft-Shelled Adlay (Coix lacryma-jobi Variety ma-yuen), a Cereal and Medicinal Crop in the Poaceae Family.

Coix lacryma-jobi, also called adlay or Job's tears, is an annual herbal plant belonging to the Poaceae family that has been cultivated as a cereal and medicinal crop in Asia. Despite its importance, however, genomic resources for better understanding this plant species at the molecular level and informing improved breeding strategies remain limited. To address this, we generated a draft genome of the C. lacryma-jobi variety ma-yuen (soft-shelled adlay) Korean cultivar, Johyun, by de novo assembly, using PacBio and Illumina sequencing data. A total of 3,362 scaffold sequences, 1.28 Gb in length, were assembled, representing 82.1% of the estimated genome size (1.56 Gb). Genome completeness was confirmed by the presence of 91.4% of the BUSCO angiosperm genes and mapping ratio of 98.3% of Illumina paired-end reads. We found that approximately 77.0% of the genome is occupied by repeat sequences, most of which are Gypsy and Copia-type retrotransposons, and evidence-based genome annotation predicts 39,574 protein-coding genes, 85.5% of which were functionally annotated. We further predict that soft-shelled adlay diverged from a common ancestor with sorghum 9.0-11.2 MYA. Transcriptome profiling revealed 3,988 genes that are differentially expressed in seeds relative to other tissues, of which 1,470 genes were strongly up-regulated in seeds and the most enriched Gene Ontology terms were assigned to carbohydrate and protein metabolism. In addition, we identified 76 storage protein genes including 18 seed-specific coixin genes and 13 candidate genes involved in biosynthesis of benzoxazinoids (BXs) including coixol, a unique BX compound found in C. lacryma-jobi species. The characterization of those genes can further our understanding of unique traits of soft-shelled adlay, such as high seed protein content and medicinal compound biosynthesis. Taken together, our genome sequence data will provide a valuable resource for molecular breeding and pharmacological study of this plant species.

De novo transcriptome sequence of Senna tora provides insights into anthraquinone biosynthesis.

Senna tora is an annual herb with rich source of anthraquinones that have tremendous pharmacological properties. However, there is little mention of genetic information for this species, especially regarding the biosynthetic pathways of anthraquinones. To understand the key genes and regulatory mechanism of anthraquinone biosynthesis pathways, we performed spatial and temporal transcriptome sequencing of S. tora using short RNA sequencing (RNA-Seq) and long-read isoform sequencing (Iso-Seq) technologies, and generated two unigene sets composed of 118,635 and 39,364, respectively. A comprehensive functional annotation and classification with multiple public databases identified array of genes involved in major secondary metabolite biosynthesis pathways and important transcription factor (TF) families (MYB, MYB-related, AP2/ERF, C2C2-YABBY, and bHLH). Differential expression analysis indicated that the expression level of genes involved in anthraquinone biosynthetic pathway regulates differently depending on the degree of tissues and seeds development. Furthermore, we identified that the amount of anthraquinone compounds were greater in late seeds than early ones. In conclusion, these results provide a rich resource for understanding the anthraquinone metabolism in S. tora.

HOS15 Interacts with the Histone Deacetylase HDA9 and the Evening Complex to Epigenetically Regulate the Floral Activator GIGANTEA.

In plants, seasonal inputs such as photoperiod and temperature modulate the plant's internal genetic program to regulate the timing of the developmental transition from vegetative to reproductive growth. This regulation of the floral transition involves chromatin remodeling, including covalent modification of histones. Here, we report that HIGH EXPRESSION OF OSMOTICALLY RESPONSIVE GENE 15 (HOS15), a WD40 repeat protein, associates with a histone deacetylase complex to repress transcription of the GIGANTEA (GI)-mediated photoperiodic flowering pathway in Arabidopsis (Arabidopsis thaliana). Loss of function of HOS15 confers early flowering under long-day conditions because elevated GI expression. LUX ARRHYTHMO (LUX), a DNA binding transcription factor and component of the Evening Complex (EC), is important for the binding of HOS15 to the GI promoter. In wild type, HOS15 associates with the EC components LUX, EARLY FLOWERING 3 (ELF3), and ELF4 and the histone deacetylase HDA9 at the GI promoter, resulting in histone deacetylation and reduced GI expression. In the hos15-2 mutant, the levels of histone acetylation are elevated at the GI promoter, resulting in increased GI expression. Our data suggest that the HOS15-EC-HDA9 histone-modifying complex regulates photoperiodic flowering via the transcriptional repression of GI.

The complete mitochondrial genome sequence of Schisandra chinensis (Austrobaileyales: Schisandraceae).

Chinese magnolia vine (Schisandra chinensis) is an economically important oriental medicinal plant that belongs to the Schisandraceae family. The complete mitochondrial genome sequence of S. chinensis was 946,141 bp in length. A total of 45 genes was annotated, including 30 protein-coding genes, 12 tRNA genes, and 3 rRNA genes. A phylogenetic tree based on the mitochondrial genome demonstrated that S. chinensis was most closely related to Schisandra sphenanthera of the Schisandraceae family.

De novo transcriptome assembly of the Chinese pearl barley, adlay, by full-length isoform and short-read RNA sequencing.

Adlay (Coix lacryma-jobi) is a tropical grass that has long been used in traditional Chinese medicine and is known for its nutritional benefits. Recent studies have shown that vitamin E compounds in adlay protect against chronic diseases such as cancer and heart disease. However, the molecular basis of adlay's health benefits remains unknown. Here, we generated adlay gene sets by de novo transcriptome assembly using long-read isoform sequencing (Iso-Seq) and short-read RNA-Sequencing (RNA-Seq). The gene sets obtained from Iso-seq and RNA-seq contained 31,177 genes and 57,901 genes, respectively. We confirmed the validity of the assembled gene sets by experimentally analyzing the levels of prolamin and vitamin E biosynthesis-associated proteins in adlay plant tissues and seeds. We compared the screened adlay genes with known gene families from closely related plant species, such as rice, sorghum and maize. We also identified tissue-specific genes from the adlay leaf, root, and young and mature seed, and experimentally validated the differential expression of 12 randomly-selected genes. Our study of the adlay transcriptome will provide a valuable resource for genetic studies that can enhance adlay breeding programs in the future.

Tracking Down a New Steroid-Hydroxylating Promiscuous Cytochrome†P450: CYP154C8 from Streptomyces sp. W2233-SM.

CYP154C8 from Streptomyces sp. has been identified as a new cytochrome P450 with substrate flexibility towards different sets of steroids. In vitro treatment of these steroids with CYP154C8 revealed interesting product formation patterns with the same group of steroids. NMR study revealed the major product of corticosterone to be hydroxylated at the C21 position, whereas progesterone, androstenedione, testosterone, and 11-ketoprogesterone were exclusively hydroxylated at the 16α position. However, the 16α-hydroxylated product of progesterone was further hydroxylated to yield dihydroxylated products. 16-hydroxyprogesterone was hydroxylated at two positions to yield dihydroxylated products: 2α,16α-dihydroxyprogesterone and 6ß,16α-dihydroxyprogesterone. To the best of our knowledge, this is the first report of generation of such products through enzymatic hydroxylation by a CYP450. In view of the importance of modified steroids as pharmaceutical components, CYP154C8 has immense potential for utilization in bioproduction of hydroxylated derivative compounds to be directly employed for pharmaceutical applications.

The complete chloroplast genome sequence of Coix lacryma-jobi L. (Poaceae), a cereal and medicinal crop.

Coix lacryma-jobi is a cereal and medicinal crop belonging to the Poaceae family. This study characterized complete chloroplast genome sequence of a Korean cultivar Johyun of C. lacryma-jobi var. ma-yuen through the de novo hybrid assembly with Illumina and PacBio genomic reads. The chloroplast genome is 140,863 bp long and composed of large single copy (82,827 bp), small single copy (12,522 bp), and a pair of inverted repeats (each 22,757 bp). A total of 123 genes including 87 protein-coding genes, 32 tRNA genes, and four rRNA genes were predicted in the genome. Phylogenetic analysis confirmed a close relationship of C. lacryma-jobi with species in the Panicoideae subfamily of the Poaceae family.

Complete chloroplast genome and 45S nrDNA sequences of the medicinal plant species Glycyrrhiza glabra and Glycyrrhiza uralensis.

Glycyrrhiza uralensis and G. glabra, members of the Fabaceae, are medicinally important species that are native to Asia and Europe. Extracts from these plants are widely used as natural sweeteners because of their much greater sweetness than sucrose. In this study, the three complete chloroplast genomes and five 45S nuclear ribosomal (nr)DNA sequences of these two licorice species and an interspecific hybrid are presented. The chloroplast genomes of G. glabra, G. uralensis and G. glabra × G. uralensis were 127,895 bp, 127,716 bp and 127,939 bp, respectively. The three chloroplast genomes harbored 110 annotated genes, including 76 protein-coding genes, 30 tRNA genes and 4 rRNA genes. The 45S nrDNA sequences were either 5,947 or 5,948 bp in length. Glycyrrhiza glabra and G. glabra × G. uralensis showed two types of nrDNA, while G. uralensis contained a single type. The complete 45S nrDNA sequence unit contains 18S rRNA, ITS1, 5.8S rRNA, ITS2 and 26S rRNA. We identified simple sequence repeat and tandem repeat sequences. We also developed four reliable markers for analysis of Glycyrrhiza diversity authentication.

Comparative transcriptome analysis reveals whole-genome duplications and gene selection patterns in cultivated and wild Chrysanthemum species.

KEY MESSAGE: Comparative transcriptome analysis of wild and cultivated chrysanthemums provides valuable genomic resources and helps uncover common and divergent patterns of genome and gene evolution in these species. Plants are unique in that they employ polyploidy (or whole-genome duplication, WGD) as a key process for speciation and evolution. The Chrysanthemum genus is closely associated with hybridization and polyploidization, with Chrysanthemum species exhibiting diverse ploidy levels. The commercially important species, C. morifolium is an allohexaploid plant that is thought to have originated via the hybridization of several Chrysanthemum species, but the genomic and molecular evolutionary mechanisms remain poorly understood. In the present study, we sequenced and compared the transcriptomes of C. morifolium and the wild Korean diploid species, C. boreale. De novo transcriptome assembly revealed 11,318 genes in C. morifolium and 10,961 genes in C. boreale, whose functions were annotated by homology searches. An analysis of synonymous substitution rates (Ks) of paralogous and orthologous genes suggested that the two Chrysanthemum species commonly experienced the Asteraceae paleopolyploidization and recent genome duplication or triplication before the divergence of these species. Intriguingly, C. boreale probably underwent rapid diploidization, with a reduction in chromosome number, whereas C. morifolium maintained the original chromosome number. Analysis of the ratios of non-synonymous to synonymous nucleotide substitutions (Ka/Ks) between orthologous gene pairs indicated that 107 genes experienced positive selection, which may have been crucial for the adaptation, domestication, and speciation of Chrysanthemum.